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An asymptotic expression is obtained for the current distribution on the outside sur- 
face of an infinitely long, perfectly conducting, hollow cylindrical antenna that is fed by an 
infinitesimally narrow circumferential gap. This asymptotic expression involves two 
series. The first series is expressed in reciprocal powers of log (2|z|/jI'%ka®), where |{z| is 
the distance from the gap, log T is Euler’s constant, k is the propagation constant, and a 
is the radius of the antenna. The second series is a similar series multiplied by 1/(k|z2|). 
The first series is dominant and its first five terms yield values for the magnitude and phase 
of the current that for even moderately thick antennas (circumferences as large as }/3) 
are accurate to within about one percent in as close as 4/3 of the gap. This is shown by 
a comparison of the values of the current obtained from these terms with the numerically 
computed values of Duncan [1962]. Asymptotic expressions for the current found in the 
literature resemble the first term of this dominant series and are accurate only at rela- 
tively large distances from the gap—except for very thin antennas. 


1. Introduction 


Consider the usual model of an infinite antenna composed of a hollow circular cylinder 
with its axis lying along the z-axis of a cylindrical coordinate system with coordinates (p, ¢, 2). 
The cylinder is cut in two at z=0 and the antenna is assumed to be excited by a sinusoidal 
generator producing an oscillating electric field H,¢’* (that is independent of ¢) across the 
gap in the cylinder formed by the cut. The voltage across the gap is 


v= ih Euz 
gap 


which may be held constant while the gap width is allowed to approach zero. As is well 
known this limiting process of allowing the gap width to go to zero leads to a singularity of 
the quadrature component of the current J(z) on the outer surface of the cylinder. The 
effect of such a gap has been considered, among others, by Infeld [1947], King [1956], Wu 
and King [1959], Chen and Keller [1962], and Duncan [1962]. 

The current on the outer surface of the cylinder is given by the formula ! 


kaV (° H® (Ba)e-7"4dy 
17 [ BH® (Ba) @) 


where 
k=2rn 
\= wavelength in free space 
a=radius of antenna 
jae 
Zo = ble 

e= dielectric constant of free space 
u=permeability of free space 
B= (kh? —y?) 

Hi” (8a) =Jo(8a) —j ¥o(8a) 

Hy” (8a) =, (8a) —j Yi (6a). 


1 See formula 13 of Duncan [1962]. 
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The last two functions are Hankel functions of the second kind, the time factor having been 
taken as e’*'. The path of integration in (1) is along C in figure 1 and is thus along the real 
axis except at y=—k where it is downwardly indented and at y= where it is upwardly indented 
to avoid the branch points of Baty=+k. The angle of 8, written 28, is chosen so that Z8=0 
for —k<y<ck. 

Duncan [1962] has recently derived expressions for /(z) that permit its computation for 
very small values of z and has obtained accurate plots of the real and imaginary parts of [(z) 
as one approaches the singularity at z=0. While this approach also makes possible accurate 
values of [(z) out to ¢=k/2| of 5 or 10, the work involved increases and the accuracy decreases 


as ¢ becomes large. 
This paper develops an asymptotic series for the current 


gi) 


I wee ~& 2 
ea ia ps eerttes ae @) 
where A, and J, are constants and 
2 
P= Spee @) 


Tt is shown that the first five terms of the first series yield a value of Z(z) that is remarkably 
accurate even for rather small values of ¢=k/z| and for thick antennas. Thus, for example 
(see figs. 2 and 3), |Z(2)| is given within an accuracy of 1 percent for all antennas thinner than 
a==27ah1=0.30 for all ¢=2r|z|\-1> 1.8 (approx.). 


2. Modification of the Integral for J(z) by Contour Integration 


This integral for 7(z) given in (1) may, by reference to figure 1, be written 


kaV 7, (| HE (Ga)e-?"" 


(Oy ct Jo BHR Bay “ 


By distorting the contour C between A and B into QO’, one sees that C in (4) can be replaced 
by OC’. It can readily be shown that in the limit of e~0 and R-- only the portions of the path 
C’ lying along branch cut 1 contribute to the integral, and hence 


La k-jo > FT (ga) H® (Bae-*) asi) 
1ea="F |, sup ay ae Ga |e - 


Decreasing the angle of the argument fa of the Hankel functions by x leads to the following 
equations: 
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HH,” (Bae~?*) = —Jy(Ba) —j Yo(Ba) = —H5” (Ba) 
A? (Bae~**) =Ji (Ba) +7 Yi (8a) =H» (Ba). (6) 


These together with the substitution y=k(1—7y) and the well-known identity 
2 
I(x) ¥ (x) —Jo(x) ¥ 1 (x) =a (7) 
makes it possible to write (5) in the form 


AVe % (° e dy 
TZ Jo n'lJo(na)+Y¥o(ne)] 


c=kz|, oka, and n=Vy (y+ 2)). (9) 


I(z)= 


(8) 


where 


The change in the contour of integration thus leaves us with a more manageable integral for 
large values of ¢ in that only the small argument behavior of the Bessel functions is important. 


3. More Tractable Expression for the Asymptotic Behavior of I(z) 


We may write (8) in form 


2Ve—% 
I(z)= 7Z, De (10) 
where 
eed a e7 “du 
= - 11 
EJ Wena) +Y ana) | Cn 
and 
| ag eee nee 
1 ey u(u+29¢). (12) 


Within the region |ya|<b, where 6 is a small positive number, the Bessel functions may be 
replaced by their small argument approximations. This amounts to setting 


Dec 1 ____1+«:(na) 
De) Sine) +Yi(na) | +(2 oe my 
T 2 


(13) 


and neglecting ¢(na). This leads to a maximum relative error in the approximation that is a 
function of 6. For 6=0.5 this maximum error is about 1 percent. 
Letting wu be the value of u for which 


| na sey ule +290] =) (14) 


we have 
UW=CF (15) 


o={ 2[ (148 1] } (16) 


It is therefore seen that for any choice of a and 6 one may make w as large as one chooses by 
taking ¢ large enough. A tabulation of C for }=0.5 and various values of @ is given in table 1. 


where 


TABLE 1 
a C for b=0.5 T rl 0.066771 Us 
0.0255 24. 96 6.34 10-8¢-1 1576¢ 104¢ 25 
0. 08_..-...- 6. 09 1.0210-2¢-1 98. 5¢ 6. 50% 2 
0. 80..._-.__- 1.19 0. 1438¢-1 7. 00¢ 0. 462¢ 0. 4625 
0. 60._-.-... 0. 34 571¢-1 1. 75¢ 116¢ 116¢ 
TW 20sc 20242 09 2, 2857-1 0. 44° 029¢ 029¢ 
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From (11) and (13) 


= sale Bes rey PS) aoa 


Since it is possible to choose an A such that |D(ya)|>(Aa|n|)“ for all |ya|>b, we may make 
use of (12) to obtain the inequality 


2Aa “o 


du 
£ ie, |n| a . ae ( ) 


al uy 7 "Dna 


which, using (15), can be written 


Fela #Doal Vor a 


Thus with the neglect of e in (17) and exponential term in ¢ we may write 


We Goes e~“du 
1M fie FY] 


where M is wu; or any quantity very large compared tol. Since the integrand of the first integral 
in (17) has nearly a constant angle, the percentage error in dropping the integral involving 


€ is certainly not much more than the maximum error in dropping ¢« in (18). 
Expansion in (20) of [(w/¢)-+27]"! by the binomial theorem leads to a series in reciprocal 


powers of ¢ 


(20) 


Tz ZyT, 21 
a (F)T., (21) 
where 
AM n-l,-u 
ge =, 7 ae 2 (22) 


At this point to insure uniform convergence of the binomial series, M is to be taken slightly 
smaller than wz, where w is the smaller of u, and 2¢. 
From (12) 


1 
= FZ no (23) 


x 
m= a (24) 


Consider now a function f of »? in the range in which uc<¢. By a Taylor expansion 


where 


orf} G2 (25) 


Since, however, 


df _dudf __ jsaf 


dnt ane du 2du 
(25) may be written 


To=soo—e Fe. (26) 
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Applying (26) to (22), we have 


un -1 e7*du 


M 
~ ntl, a 
a ° 1+ (Eog Ey eS i 35 Ting EY Ey | 
T 


The application of (26) in the outer part of the range 0<uw<M will not be very accurate, but 
since the integrand drops off as e~“, one may ignore this. Integrating the last integral by parts 
and neglecting the term containing e~™ in the integrated part, we have 


Tye Tt L(t) TL i el (27) 


where the superscript zero indicates that 7 is to be replaced by 7. 
If in (21) one neglects terms of (1|¢)? and higher, then one may write that 
Tx Totgeh 
By (27) this becomes 
T= T3+3-[20%—T3), (28) 


Equations (10) and (28) provide us with an asymptotic expression for J(z) involving 
integrals of the form 


es unt a 
? me log 1+(Hog !nie* ; 
we yt -1 see 
=. ee y 
1+-(Zlog UT (29) 
where use has been made of (24) and where 
2,2 
=< . (30) 
4. Asymptotic Series for I(z) 
Integrating (29) by parts using the function 
A(2) =5tarctan x (81) 
we have 
bens 1 \y, ¢¥ a 1 
a Tr=| ure*A (Flog uri) | +f (u™—nu"")e~"A € log uri) du. 
Since A(e log urj) 00 as u 0 and M is large, this reduces to 
T° —W,—nW,-1-+Mre-™A é log Mri) 
2W,—mW n-1 (32) 
where 
tt 1 : 
Waar ure-"A (= log uri) du. (88) 
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Collecting our results thus far, we have from (10), (28), and (82) that 


9 -i# } 
ot _ iz LW, 2W.—W.) (34) 


Corresponding to the expansion of the arctangent in reciprocal powers of the argument, 
one has 


(—1)*2/1 \ = Ok+1) 
A(= log urj = LI (Zlog wri) . (35) 


i) 
This series converges for 


pk 


ieee 
etd 


and thus will converge when 


—log ur>Baror U<c0.066771, (36) 


Let wu; be the smallest of the quantities u.=Cf, 2f, and 0.06677! then (35) converges 
uniformly for 0<u<M<us, and with neglect of terms of order e~” 


= Come 2k+2Q 
Wied, QL T re (37) 
with 
M Nyp—tU 
ps ure “du (38) 


0 (og uri yr 


This follows from the fact that a uniformly convergent series may be integrated term by term. 

It is here understood that M is to be taken only negligibly smaller than u;. It is seen from 

table 1 that u3, and hence M, can be taken as large as one desires by increasing ¢. 
Equation (38) may be written 


ure "du 
(log sa |, [1+ 008 exp (39) 
log 77 


Lua= 


Using Taylor’s series with a remainder term to expand the denominator, we have (see Widder 


[1961]) 
log U ey N ti —log wr? 
Ear = 25 ery) ( log =) +By (40) 


where {s} as a superscript is a pseudoexponent defined by the equation 
( (2k-+1)(2k+2) ... (2k-+s), s>0 
(2k+1)=< 1, s=0 (41) 
[(2k)(2k—1) ... (Qk+s8+1)]7}, s<i0 


and 


ae —log u_ ‘| 
fee (2k-+1) eee Te ee ee 7 dt (42) 


(8) | 2k+ N+2 


Letting ¢=—log v/log 77 in (42), one may write 


~ Nidog 77) *F log or) (log orj\*tNt2 > (48) 


log r7 


ee ee 
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So ES rn 


7 


On substituting from (40) into (89), one obtains the equation 


{s} 
I= gga | 2 EE gS + Rte) | (44) 
where 
Cu=(—1)" [" ure—“(log u)sdu (45) 
and : 
Ri i "wen" Ryrd. (46) 


It will be shown later in section 5 that as f-« and hence 7-0 the remainder term y(n) 
approaches zero as |log zj/-“*?, This insures that (44) is a legitimate asymptotic expansion 
of Luis 

It is convenient to replace (45) by the equation 


Cle (-1y[" ue“ (log u)sdu (47) 
0 


which involves again only the neglect of a term proportional to e~™ and to write (44) in the form 


: -1 & | kt Cus 
mE Togry =i st log rj) 


(48) 


to indicate an asymptotic series that must be terminated at some finite but unspecified term. 
The remainder term is then very roughly indicated by the first neglected term. A precise 
determination of the overall error is made later in section 7 by comparing the results obtained 
from our asymptotic series with those obtained numerically by Duncan [1962]. 

Integrating (47) by parts, we have 


(—1)0 a] urtte» COROT 


= (~ [u™t!—(m+1)u"e—*(log u)st du 


gaat _lu TS Jd 
or 
Crt, i= (N41) C;, o41—- (8s +D Ci, ¢- 
Thus 
Cy Gens}, s—8Cn-1, s—1 (49) 
and, in particular, 
Oe Cos—sCo, s-1 (50) 
Cr5=20C1.—8C), s—1 
=2C,—38C, s—1 + 8(s—1) Cy, s~2s (51) 


From (37) and (48), since in (87) W, can be approximated by a truncated series, 


Wotse [4W,—2W,—W)] 


(0) (—1)*"} aa 2k+1 (co) (2k+1)'s) oat 
k=6 2k+1 log =) s=0 s! (log Toe al oT 2¢ 


and thus, finally, by (84), (41), (50), and (51) 


= 


I (40,,—20.— 0) | (52) 


ae If (@) (@) (— 1)*+3(2he4-2) 8-72 ee 
He) SS aoe RT [OR COatee—-N0-) | 3) 
where 
a on nn {, ” e- “(log u)tdu. (54) 
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Equation (53) may be written 


= x B 
ig : Terry | 
Oe > ea fee, pe Hy (log pg)” ee) 
where A’s and B’s are constants and —log 77 is replaced by log pr; thus 
—2) 
P=Tar (56) 


These are the same equations as (2) and (3) in the introduction. Vainshtein [1959] has obtained 
a similar type of expansion for a related integral. 

The second series in (55) is not, strictly speaking, a part of an asymptotic expansion of 
I(z) since terms in 1/¢ fall off faster than any power of (log pf)~*. It can, however, be used 
to increase the range of ¢ over which (55) is a valid approximation. 

By direct comparison of (53) and (55) we see that 


A,=C 
A,=—0O; B.=} Ch 
2 
A=—5 Q+C; By=—(Q+C) (57) 
1 3 
A,=7?0,—C; Bay Corsi ts C; 


As= Oy 20 Ort Cy By=202(Cy-+-C,)—6C—2C, 


5. Magnitude of the Remainder Term Rx(n) 
If uw is in the range 0<u<77™, then from (43) 
2E+N +2 
log r+ 54 


\ v N+1 44 
[S2-|k 


Since [log v7j| in (48) for real v is always greater than 7/2, we have 


(Qh-- 1) NFU Q2R+N42 Hog y|NH 


(Qh--1) (Nt ; 
(NEI)! logrgtt «(58) 


[Rwl< WH 
Nilog I" Toe 483 


ip eteer) (N+ log 79/** log uit (59) 


2k+N+2 
(w+n!(§) 
even outside the above range. 
These inequalities applied to (46) yield an upper bound for Ry(n) expressed by 
kei) tu Qek+N+2 


, (2 iti ee : 
[Rv@)|<~ aay Tog RF Jo ure~“log ul¥ du 


o 7 4[2ktl 
4 Hlog flog 7g|** hee ure-*| log udu b (60) 


TON 2k-FN-F2 
® 
— 1/2 


The second integral approaches zero as e~* ” and thus for small 7 (large ¢) is negligible. One 
sees therefore that Ry(n) goes to zero as 
constant (61) 


flog 79°F" ” 
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6. Computation of the C,, 
From. (54) 


C,,=(—1)”" i, e~ “(log u)"du=P,,+(—1)"Q, 
where ° 
P,=(—1)" [ e-“(log u)"du>0 
and : 
o= { ” e-4(log u)"du>>d- 
In (68) let z= —log u; then : 
Pra [ere caserie 
=, (—1)? 


® a eae 
pe ne (atid, —n! eee Ae 
a p! x os Zplp ty 


and one obtains a convenient alternating series for P,. 
To obtain the Q, we set 


MW 
Q,= |) e "dog u)"dute, (Md) 
Jl 


(62) 


(63) 


(64) 


(65) 


(66) 


where the truncation error e, (14) in using VM rather than infinity as the upper limit is given by 


€,(M) =|" e“(log u)"du=e sad Me a [ og +41 og ( +a) | 


Now log +5) S47 <—~_ for all x=>0 and hence e,(M) is less than 


E,(M) = 


M 
aa e*(M log M+-2)""x 


which on successive integration by parts reduces to 


5 Atlee MM) 


E,(M)= I 


= (log Af)"e-* 


n n(m—1) nr! 
E 'MlogM* (Mlog Me" °° * * (flog | 


Setting n=5 and\M=16, we have 
£;(16)=(2.77)8e-"[1.124]=2.08 X 1078, 
For smaller n the error is even less. Thus, we may with negligible error set 
16 
Qu=f e “(log u)"du, n=0,1,2,... 
and evaluate these integrals by numerical quadrature. Actually, one knows that 


C= | "e-vdu=1 
0 


and from standard integral tables that 
a={ e~* log udu=y7=Euler’s constant. 
0 


Hence numerical integration is not required in (69) for n=0 and n=1. 
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(67) 


(68) 


(69) 


(70) 


(71) 


A table of the C,, and the A, and B, obtained from them by (57) is given below: 


TABLE 2 
n | C, | An RB, 
1.0000 
0. 5772 1. 0000 
1.9781 | —0.5772 0. 5000 
5.4449 | —1.3118 —1. 5772 
23. 561 0. 2520 —0. 2360 
117. 84 3. 9969 8. 3743 


7. Accuracy of Results Obtained Using Asymptotic Series 


Since several approximations were made in deriving the asymptotic series (55) for the 
current on an infinite cylindrical antenna, it is important to check the results against the very 


accurate results obtained by Duncan [1962] by numerical integration. 


For this purpose only 


the first five terms of the A coefficient series were used. The computations were made for rather 
thick antennas, a=0.02, 0.08, 0.30, 0.60, and 1.20, and some of the results are shown in 


figures 2 to 10. 
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Figure 2. Magnitude of the current for a=0.08. 
“Asymptotic” refers to first five terms of the A,scries in (55). 
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Figure 3. Magnitude of the current for e=0.30. 
“Asymptotic” refers to first five terms of the A series in (85). 
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Fiaure 4. Magnitude of the current for e=0.60. 
“Asymptotic” refers to first five terms of the A series_in (55). 
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Ficure 5. Magnitude of the current for a= 1.20. 
““Asymptotic”’ refers to first five terms of the A series in (55). 


The results for a=0.02 are not given since at a=0.08 (see figs. 2, 6, and 9) the results are 
already quite accurate even as close as a quarter wavelength from the gap. For thinner 
antennas the agreement would be even better. One observes (see figs. 4, 7, and 10) that for 
a=0.60 the asymptotic series begins to give accurate results only when ¢ is greater than 3.5 
or 4. For a circumference approaching one wavelength (a=1) the asymptotic expansion 
ceases to be even approximately valid for determining the magnitude of the current as can be 
seen from figure 5 (a=1.2). Surprisingly enough figure 8 shows that the phase angle of the 
current for a=1.2 is somewhat better predicted, being less than 20 deg off from the true value. 
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Fieure 6. Phase angle of the current for e=0.08. 
“‘Asymptotic”’ refers to first five terms of the A series in (55). 
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Fieurs 7. Phase angle of the current for a=0.60. 
“Asymptotic” refers to first five terms of the A series in (55). 
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Figure 8. Phase angle of the current for a=1.20. 
‘“Asymptotic” refers to first five terms of the A series in (55). 
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Ficure 9. Imaginary part of the current for a=0.08. 
“ Asymptotic” refers to first five terms of the A series in (55). 


429 


5.0 LEGEND 


® EXACT VALUES 


4.0 
a =0.60 i 
* \ 
3.0 / \ = FIRST TERM 
a 
bs I, 
Q oy EXACT e 
2.0 \ \ > 
wo 
2) n 
[o) 
\ Ose 
le m 
Ag : 3 
rls \ 2 
x 0 — o & 
bE TO0~,__8.0\.4— 9.0 10.0 0 ————_-_» 
g C 
a @ 
B-1.0 ‘\ 
-05 
i) 
-2.0 
Wes 
\ -1.0 
% 5.0 
-3.0 / ~~ 


ee. 


Figure 10. Real part of the current for a=0.60. 
‘Asymptotic’ refers to first five terms of the A series in (55), 


The first term of the asymptotic expansion (55) taken alone yields the asymptotic formula 


y) V —jklzl 
0 log J] val 
Ika? 


which is accurate for antennas of the thickness here considered only for rather large ¢=k{z). 
This is shown in figures 9 and 10. Even for the a=0.08 case the plot of the first term fails to 
give close agreement even at a ¢ of the order of 10. These two figures illustrate, in addition, 
the degree to which the first five terms give the real and imaginary parts of the current. Since 
the accuracy is approximately the same for each part, we have given only the imaginary part 
for a-=0.08 and only the real part for a=0.60. In the latter, of course, the agreement is not 
too good. 

Formulas similar to (72) have been obtained by others. Chen and Keller [1962] obtain the 
same formula but with the argument of the log the negative of that in (72). Northover [1958] 
derives a formula in which the argument of the log is 2/2ka?. Hallén’s [1956] formula agrees 
with (72) except that T is to the first power. All of these formulas are correct in that they lead 
to the right functional behavior for s->. It has been checked, on the other hand, that in 
the case a=0.08, at least, formula (72) is more accurate for finite ¢ than any of these variations 
of it. 

We have also determined that for a=0.08 and 0.60 the error remaining after taking the 
first five terms of the A coefficient series is given for intermediate ¢ very nearly by the first 
two or three terms of the B coefficient series of (55). 
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